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ABSTRACT 

Recently the initial supersonic relative velocity between the dark matter and baryons 
was shown to have an important effect on galaxy formation at high redshift. We study 
the impact of this relative motion on the distribution of the star-forming halos and on 
the formation redshift of the very first star. We include a new aspect of the relative 
velocity effect found in recent simulations by fitting their results to obtain the spatially- 
varying minimum halo mass needed for molecular cooling. Thus, the relative velocities 
have three separate effects: suppression of the halo abundance, suppression of the gas 
content within each halo, and boosting of the minimum cooling mass. We show that 
the two suppressions (of gas content and of halo abundance) are the primary effects 
on the small minihalos that cannot form stars, while the cooling mass boost combines 
with the abundance suppression to produce order unity fluctuations in stellar density. 
We quantify the large-scale inhomogeneity of galaxies, finding that 68% of the star 
formation (averaged on a 3 Mpc scale) is confined to 35% of the volume at z = 20 
(and just 18% at z = 40). In addition, we estimate the redshift of the first star to be 
z ~ 65, which includes a delay of Az ~ 5 due to the relative velocity. 

Key words: dark ages, first stars, large scale structure of the Universe 



1 INTRODUCTION 

In the present era of "precision cosmology" and rapidly 
advancing observational capabilities it is important to 
make precise theoretical predictions for future observations. 
Among the major goals of observational cosmology in the 
near future are to collect data on structure at high redshifts 
(including the first galaxies), detect the 21-cm line of inter- 
galactic hydrogen, and study the cosmic reionization history. 
A deep understanding of structure formation on small scales 
and at high redshifts is crucial for making reliable predic- 
tions that will help us explore this observational frontier. 

The linear perturbation theory of structure formation 
in the framework of the flat ACDM model is well under- 
stood. It allows us to follow the evolution of structure start- 
ing from tiny perturbations. The large-scale perturbations 
are 0(1O -5 ) of the b ackground quantities at cosmic recom- 
bination, z ~ 1100 ()Komatsu et al.ll2010l h and may have 
been produced during an early period of inflation. Structure 
on the smaller scales on which halos form evolves nonlin- 
early. In order to make reliable predictions, it is important 
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to verify when we can trust the results of the linear pertur- 
bation theory and on which scales the nonlinear effects must 
be accounted for. 

Linear theory separates different scales, so that each 
density perturbation mode at a given wavenumber k evolves 
independently. Thus, nonlinear terms that couple the large- 
scale velocity to the small-scale density perturbations are 
neglected in linear perturbation theory. Ho wever, recently it 
was shown jTseliakhovich fe Hiratal l20TCh that such terms 
might be of the same order of magnitude as the linear terms 
exactly at the time and on the scales on which the first bary- 
onic objects formed. Specifically, the photon-baryon cou- 
pling before recombination left the dark matter and bary- 
onic fluids with large relative velocities. These velocities im- 
pede the gravitational perturbation growth on small scales, 
leading to a sp atially-variable suppression in the abun- 
dance of halos (1 Tseliakhovich fc Hiratal l201Ch . Moreover, 
halos that later form cannot ac crete the gas as it shoots 
past the collapsing dark matter (Dala i. Pen fc Seliakl |20ld : 
iTseliakhovich! Barkana fc Hiratall2010l h 

In this paper we study the impact of the relative ve- 
locities on the distribution of the star-forming halos at high 
redshift and on the redshift of formation of the very first star. 
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In particular, we include an aspect of the relative velocity 
effect that has not been previously accounted for, and which 
is critical for understanding the overall impact of the veloci- 
ties on the distribution of star formation. Recent small-scale 
numerical simulations (|Stacv et al.ll20icil ; IGreif et al.ll201lT ) 
found that the relative velocity substantially increases the 
minimum halo mass in which stars can form from gas that 
cools via molecular hydrogen coolin g (The effec t of th e ve- 
locities has also b een simulated by iMaio et alj (|201lT ) and 
iNaoz et all (|201ll )). 

This paper is organize d as follows. In Section 2 we 
brie fly review the results of ITseliakhovich fc Hiratal (|2010l 1 
and ITseliakhovich. Barkana fc Hiratal l|201(t ). In Section 3 
we summarize the results of recent simulations that include 
the nonlinear effect of the relative velocity on the formation 
of the first stars via molecular cooling. We use the simu- 
lation results to find the behavior of the minimal cooling 
mass versus redshift and magnitude of the relative velocity. 
In Section 4 we study in detail the probability distribution 
of the gas fraction in halos at high redshift, separating out 
and comparing the importance of the various effects of the 
bulk velocity. In Section 5 we then estimate the redshift of 
the very first star accounting for the relative velocity effect. 
Finally, in Section 6 we summarize our results and also give 
a comp lete list of differences compar e d to three previous 
paper s: ITseliakhovich fc Hiratal (fcOloMDalaL Pen fc Seliakl 
i|2010l ). and ITseliakhovich. Barkana fc Hiratal i^Old ). 

Our calculations are carried out in a flat ACDM uni- 
verse with cosmological parameters taken from the 7-year 
WMAP results (WMAP7+BA O+J/q maximum likelihood 
fit from iKomatsu et al.l l|2010l )): the dark matter density 
today Sl Ci o = 0.2265, the baryon density Qb,o = 0.0455, 
the vacuum energy density J1a = 0.728, the Hubble con- 
stant Ho — 70.4 km s _1 Mpc -1 , and the spectral index 
n s = 0.967. We normaliz e the power spectru m to give a 
present value of as = 0.81 l|Komatsu et alj201Ch. We use the 
CAM B-sources linear perturbation code (|Lewis fc Challinorl 
120071 ) to generate initial conditions at recombination (specif- 
ically, at z = 1020 and z = 970 in order to obtain the needed 
derivatives) . 



2 REVIEW OF THE RELATIVE VELOCITY 
EFFECT 

In this section we briefly review the non-linear effect of 
the relative velocities between the baryons an d dark mat- 
ter, as discussed in |Tseliakhovich fc Hiratal |2010) and 
ITseliakhovich. Barkana fc Hiratal (|2010h . the latter of which 
we closely follow in our subsequent calculations. 

The initial conditions at recombination include signif- 
icant relative velocities between the baryons and the cold 
dark matter (which we denote «b c )- Before the baryons kine- 
matically decouple from the radiation (around z = 1100), 
they are carried along with the photons, while the dark 
matter moves according to the gravitational growth of fluc- 
tuations which has been advancing since matter-radiation 
equality (2 ~ 3200). At decoupling, the baryonic speed of 
sound drops precipitously, and the relative velocity then be- 
comes a substantial effect. 

In the standard picture of Gaussian initial conditions 
(e.g., from a period of inflation), the density and the com- 



ponents of relative velocity are Gaussian random variables. 
The velocity and density are spatially correlated (at differ- 
ent points) since the continuity equation relates the veloc- 
ity divergence to the density. Indeed, this equation gives an 
extra factor of 1/k in the velocity (where k is the wavenum- 
ber), making the velocity field coherent on larger scales than 
the density. Specifically, velocity fluctuations have signifi- 
cant power over the range k ~ 0.01 — 0.5 Mpc -1 . 

The relative velocity is thus coherent on scales smaller 
than ~ 3 comoving Mpc. We therefore analyze probabil- 
ity distributions in such coherent patches, and refer to the 
uniform relative velocity within each patch as the "bulk" 
or "streaming" velocity. The magnitude of the bulk veloc- 
ity in each coherence patch at recombination is distributed 
according to a Maxwell-Boltzmann distribution function: 



Pv bc (Vbc) = 



M- 



3/2 



47ru bc exp 



'2<j2 



(1) 



where a Vbc ~ 30 km sec -1 is the root-mean-square velocity 
at recombination. Just like any peculiar velocity, the bulk 
velocity Ub c decays as (1 + z) with the expansion of the 
universe. In addition to the bulk velocity, within each patch 
there are small-scale peculiar velocities of the baryons and 
dark matter related to the evolution of perturbations (and 
formation of halos) within the patch. 

As was shown in the above references, inside each co- 
herent region the linear evolution equations for density and 
velocity perturbations are modified. For example, on small 
scales the nonlinear term in the continuity equation that 
couples the local density to the velocity field, a~ 1 v ■ V<5, 
is comparable to linear terms such as the velocity term 
a _1 V ■ v. The leading contribution of the nonlinear term 
comes from the bulk motion (a _1 Vb c • V<5) and this contri- 
bution is then linear in terms of the perturbations within the 
patch. As a result, the evolution equations for the perturba- 
tions inside a coherent patch are still linear but dependent 
on the bulk Vh c - The resulting velocity-dependent terms were 
previously neglected but must be included when structure 
on small scales and at high redshifts is considered. 

The relative velocity effect is particularly important for 
the formation of the first stars and galaxies. As the first 
baryonic objects try to form, they must do so in a moving 
background of the dark matter potential wells. This rela- 
tive motion means that the dark matter's gravity must work 
harder in order to trap the baryons. As a result, the forma- 
tion of the first bounded baryonic objects is delayed. The 
effect, though, is less relevant for structure formation to- 
day, since the relative velocity decays with time while the 
typical mass of galactic host halos increases. However, the 
relative motion may shift slightly the positions of the BAO 
peaks an d produce a unique signat ure in the bispectrum of 
galaxies l|Yoo. Dalai fc SeliaJdl2oTll ). 



3 CALIBRATION OF THE MINIMUM HALO 
COOLING MASS WITH SIMULATIONS 

The formation of the first baryonic objects (in particular the 
first stars) was an important milestone in the history of the 
Universe. It marked the transition between the cold, neu- 
tral, metal-free universe (the epoch called the "dark " cos- 
mological ages that started right after recombination) and 



© 0000 RAS, MNRAS 000, 000-000 



Impact of Relative Motion on the First Stars 3 



the modern ionized, hot, and metal-rich universe. The for- 
mation of the very first stars is expected to be relatively 
simple; this is due to the primordial chemistry before stars 
produced heavy elements, and the simplified gas dynam- 
ics in the absence of dynamically-rel evant magnetic fields 
and feedback from lum inous objects (|Tegmark et al.lll997l ; 
iBarkana fc Loebll200ll ). 

Since molecular hydrogen line emission is the lowest- 
temperature coolant in metal-free gas, the first stars are ex- 
pecte d to have formed in h alos with total mass above ~ 10 5 
M (|Tegmark et alJll997h . More generally, if the mass of 
a dark matter halo is higher than a threshold referred to 
as the minimum cooling mass (M coo i), the collapsing gas is 
heated to a high enough temperature that it emits radia- 
tion. It then cools and condenses, allowing a star to form. 
The threshold can also be described as a minimum circular 
velocity (V coo i) via the standard relation V c = \J GM/R for 
a halo of mass M and virial radius R. 



This scenario of the earliest star formation has been 
confirmed by numerical simulations using both Adap- 
tive Mesh Refinement (AMR) and Smooth Particle Hy- 



drodynamics (SPH) codes (e.j 


I.. iFuller & Couchman 2000; 


Abel. Brvan & Norman 


2002 


; iBromm. Co 


ovie & Larson 


2002 




Yoshida et al. 2003 


; iReed et all 12005; 


Yoshida et al. 


2006 


Turk et al.l 201lf). All these simulations did not ac- 



count for the initial relative velocities between the baryons 
and the dar k matter. We now summarize two recent SPH 
simulations jStacv et al.ll2010l ; IGreif et al.ll201ll ) that stud- 
ied the impact of the relative streaming velocity Vbc on the 
formation of the first stars. 



Numerical simulations face a great difficulty at high red- 
shift, since they must resolve the then-typical tiny galaxies 
while at the same time capture the global galaxy distribution 
which is cha racterized by strong fluc tuations on surprisingly 
large scales |Barkana fc Loebll2004l ). The relative velocities 
are correlated up to scales above 100 Mpc, and they are im- 
portant at high redshifts where star formation is dominated 
by very small halos. Cosmological simulations that cover this 
range of scales are not currently feasible. 

However, numerical simulations are the best tool for 
studying the complex, non-linear formation of halos on small 
scales. The scales relevant to the formation of the small ha- 
los that host the first stars are well below the coherence scale 
of the relative velocity field. Therefore it is possible to sim- 
ulate halo formation in small patches of uniform Vbc- The 
simulations yield the mass reached by a halo when it first 
allows a star to form, i.e., when it first contains a cooling, 
rapidly-collapsing gas core. The results show a substantially 
increased halo mass in regions with a significant relative ve- 
locity. This is a different effect from the suppression of the 
amount of gas, which implies a smaller number of stars in 
the halo at a given time; instead in this case there is a sub- 
stantial delay in the formation of the first star within the 
halo. Moreover, this effect is not simply related to the total 
amount of accreted gas, since in the cases with a bulk veloc- 
ity, even if we wait for the halo to accrete the same total gas 
mass as its no-velocity counterpart, it still does not form a 
star (even within the now deeper potential of a more mas- 
sive host halo); the delay is substantially longer than would 
be expected based on a fixed total mass of accreted gas. 
Instead, it appears that the explanation lies with the inter- 
nal density and temperature profiles of the gas, which are 



strongly affected by the presence of the streaming motion. 
A plausible explanation for the resulting delay in star for- 
mation is that the first star forms from the gas that would 
have accreted early and formed the dense central cores in 
which stars form; this gas tries to accrete early (when Vbc 
is still very large) into a still-small halo progenitor, so it is 
affected most strongly by the suppression of gas accretion 
due to the bulk velocity. 

The simulations yield a minimum halo cooling mass at 
various redshifts, so we fit the results to find the dependence 
of the minimum halo mass on the redshift of collapse and on 
the bulk velocity, v^ c , in the patch. This will then allow us 
to study the effect of the relative velocity on the formation 
of the first stars using statistical methods that average over 
larg e cosmologic al regions tha t cannot b e direc tly simulated. 

IStacv et all (l2010t ) and IGreif et all l|201ll) state appar- 
ently contradictory conclusions, one claiming a negligible 
effect on star-forming halos and the other a large effect. In 
order to meaningfully compare their results, it is important 
to put them both on the same scale. We express the cooling 
threshold as a halo circular velocity, since simulations (cited 
above) without the bulk velocity find an approximately 
redshift-independent threshold V C ooi,o; this is naturally ex- 
pected since molecular cooling turns on essentially at a fixed 
gas temperature, and the halo circular velocity determines 
the virial temperature to which the gas is heated. Thus, the 
limit of zero bulk velocity simply gives a fixed threshold 
Vcooi.o- When we add the relative velocities, in principle the 
minimum circular velocity in a patch may be a separate func- 
tion of two parameters, the redshift z and the bulk velocity 
at halo formation Vbc(z). The history of Vbc at earlier red- 
shifts cannot introduce additional parameters, since given 
both z and Vbc(z), the full history of Vbc is determined, i.e., 
at any other redshift z' , Ubc(z') = Vbc{z) x (1 + z')/(l + z). 

Consider now the limit of a very high bulk velocity, 
Vhc{z) ^> Vcoo^o, so that the effect of Vcooi, is negligible. 
For simplicity, consider for a moment a constant Vhc versus 
redshift, fixed at its final value Ubc(-z) at the halo formation 
redshift 2. In this case there is only one velocity scale in the 
problem. As in a Jeans mass analysis, in the reference frame 
of a collapsing dark matter halo with a circular velocity V c , 
clearly gravity will be able to pull in the gas (which streams 
by at the velocity Vbc{z)) if V c > Vbc(z). Now, in the real case 
where Vbc(z') is higher during the formation of the halo, we 
would expect to get a threshold that is higher than Vbc(z), 
but by a fixed factor, because the physics is scale-free: on 
one side, Vbc scales in a simple way with redshift, and on the 
other side, halo formation (in the high-redshift, Einstein de- 
Sitter universe) also scales in a simple way, as we know from 
spherical collapse; e.g., turnaround for a halo that forms at 
redshift z always occurs at z where 1 + z = 1.59(1 + z) so 
that Vbc{z') = 1.59ubc(z). The only new scale that enters is 
from «bc at recombination, but as long as we consider halos 
that form long after recombination, this should be insignifi- 
cant. 

Thus, the threshold circular velocity Vcooi should change 
continuously between two limits, Vcooi = Vcooi, when 
«bc(z) < Vcooi, 0, and Vcooi = av hc (z) when v hc (z) > V C ooi,o 
(in terms of a fixed, dimensionless parameter a). When Vcooi 
is expressed as a function of Vh c {z), there is no additional 
dependence on z in these two limits, so we might naturally 
expect this to be true in the intermediate region as well. 
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Indeed, the above argument suggests more generally that 
halo formation and Vbc(z) scale together so that the effect of 
the bulk velocity should not depend separately on redshift; 
also the effect of molecular cooling is a redshift-independent 
threshold. Thus, when both effects act together, the result 
should still depend on just one parameter. 

We expect the dependence on velocity to be smooth and 
well-behaved for vector Vbc(z) near zero, i.e., as a function 
of the velocity components. This suggests a quadratic de- 
pendence on [vbcf^)] 2 = [vbc(-z)] 2 rather than, e.g., a linear 
dependence on Vbc(z)- We thus propose a simple ansatz for 
the minimum cooling threshold of halos that form at redshift 
z: 

V cool (z) = {Kooi.o + [m; bc (,3)] 2 } 1/2 . (2) 

The dependence of the circular velocity V coo \ on redshift only 
through the final value i>bc(z) implies that the star- formation 
threshold in a patch with a statistically rare, high value of 
Vbc at low redshift is the same as the threshold in a patch 
with the same (but now statistically more typical) value of 
Vbc at high redshift. This should be the case during the era 
of primordial star formation, before metal enrichment and 
other feedbacks complicate matters. 

We summarize the results of the two simulations to- 
gether with the best fits to each of them (with V coo \fi and 
a as free parameters) in Figure [l ] (top panel) . We obtain 
four data points from lStacv et alTlj2010l ) with non-zero ve- 
lociti es (and two more at Vb c (z) — 0), and three points 
from iGreif et all ( |201lf > (plus three more at Vbc(z) = 0). 
The best-fit parameters are: (1) V cooi.o = 3.64 km sec~ 
and a = 3.176 for the results of Istacv etakl i201Ch: (2) 
V cool .o = 3.786 km sec -1 and a = 4.707 for IGreif et ail 
l|201ll ). 

We note that despite the small numbers of halos, we 
would not necessarily expect as large a scatter in the mea- 
sured Vcooi(^) as in other measurements of halo properties; 
for example, in a sample with a large number of halos of 
various masses at each redshift, we would expect a large 
range of redshifts for the first star formation within a halo, 
but if we only take halos that first formed a star at a given 
redshift z, their masses at z might span a narrow range, all 
near the minimum cooling mass for that redshift (since any 
halo well above the cooling mass at z would already have 
formed a star earlier). In any case, our ansatz fits each set 
of simulation results reasonably well, but there is some scat- 
ter a nd also a sy s tema tic difference between the two sets 
(with IGreif et alj (|201ll} indicating a stronger effect of the 
bulk velocity). Due to the small number of simulated halos, 
it is difficult to separate the possible effects of different nu- 
merical resolutions, other differences in the gravitational or 
hydrodynamical solvers, and real cosmic scatter among ha- 
los. Given the systematic offset, we do not simultaneously fit 
both sets of points, but instead average the best-fit parame- 
ters of the two SPH simulation sets. We mostly use this fit, 
which we refer to as our optimal fit, in the following sections: 

V cool (z) = {(3.714 km/s) 2 + [4.015 • v bc (z)] 2 } 1/2 . (3) 

There is some discrepancy in the value of V CO oi,o found in 
AMR and SPH simulations. In order to test the full current 
uncertainty range including different types of simulations, 
we also consider the avera ge value V^ooi.o ~ 4.2 km sec -1 
found in AMR simulations (|Yoshida et al.ll2006l : iTurk et all 
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Figure 1. Top panel: The minimum halo circular velocity for 
gas cooling via molecular hydrogen versus the bulk velocity v^c (z ) 
when th e halo virial i zes. D ata are taken from lStacv et all j2010h 
(•) and IGreif et al ] Hon]) (□). We show our fits to each set of 
simulation results (dot-dashed and dashed, respectively). We also 
show our "optimal" fit to SPH simulations (thick solid line), the 
"fit" to AMR simulations (regular solid line), and the case of no 
streaming velocity (dotted line, based on our optimal fit). The 
vertical solid line marks the root-mean-square value of i>bc( z ) at 
z = 20. Bottom panel: We show the minimum halo mass for 
molecular cooling versus redshift, in a patch with the root-mean- 
squarc value of Vb c (z) at each redshift z, for each of the fits from 
the top panel; in particular, we show (dotted line) the case of no 
relative motion based on our optimal fit (i.e., V coo \ = V coo i,o = 
3.714 km sec" 1 ). 



l201ll ). Thus, we combine this value of V^ool with a from 
our optimal fit to obtain what we refer to as a "fit" to AMR 
simulations. In other words, we assume that the discrepancy 
between the two simulation methods is only in the cool- 
ing process (due to systematic entropy differences in dense 
cores), but that they would agree on the effect of the bulk 
motion. Regardless of which fit we use, Figure [1] shows that 
the relative motion has a large effect on the minimum cir- 
cular velocity. 

The implications for the minimum cooling mass as a 
function of redshift are also shown in Figure [T] (bottom 
panel). In a patch with no relative motion, the mass drops 
rapidly with redshift, since at higher redshift the gas density 
is higher and a given halo mass heats the infalling gas to a 
higher virial temperature. However, in a region at the root- 
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mean-square value of v^ c the higher bulk velocity at high 
redshift implies that a higher halo mass is needed for effi- 
cient molecular cooling. In particular, at redshift 20 a patch 
with Vhc = will form stars in 3.6 x 10 5 Mq halos, while 
a patch with the root-mean-square value of Vhc has a mini- 
mum cooling mass of 6.0 x 10 Mq according to the optimal 
fit, or a range of (4.8 - 7.3) x 10 s Mq from the other fits. At 
z = 60 these numbers become 7.2 x 10 4 M Q , 7.0 x 10 s M , 
and (4.1 — 10.3) x 10 5 Mq, respectively. In patches with low 
bulk velocity we expect stars to form earlier, since the halos 
with lower masses are more abundant and form earlier in 
the hierarchical picture of structure formation. This is the 
basis of the discussion that follows. 



4 GAS FRACTION IN THE FIRST BOUND 
BARYONIC OBJECTS 

4.1 Global average 

The population of gas-filled halos at high redshift divides 
naturally into two major categories. The first category con- 
sists of large halos in which the gas can cool (via molecu- 
lar hydrogen cooling); these are presumed to be the sites 
of formation of the first stars, and are obviously most im- 
portant since the stellar radiation is in principle observable, 
and it also produces feedback on the intergalactic medium 
and on other nearby sites of star formation. Also interest- 
ing, though, is the second category, namely the smaller ha- 
los ( "minihalos" ) in which the gas accumulates to roughly 
virial density and yet cannot cool. The latter may affect the 
epoc h of reionization by acting as a sink for ionizing photons 
fe.g..lHaiman. Abel fe Madaull200ll; iBarkana fe Loebl|2002i ; 
llliev. Scannapieco fe Shapiro! 12005^ ICiardi et all I2005T I and 
may gen erate a 21-cm signal from collisional excitation of 



Hi (e.g.. llliev et"alll2003l ; iFurlanetto fe ObJl2006l V 

In this subsection we apply the result we found for the 
minimum cooling mass to find the redshift evolution of the 
gas fraction in these two categories. In the following subsec- 
tions we explore the probability distribution function (PDF) 
of the gas fraction, beginning with its dependence on the 
bulk velocity. In addition, though, in each patch of coherent 
velocity the mean density is slightly different, varying as a 
result of random density fluctuations on scales larger than 
the patch size. We thus also study the full PDF as deter- 
mined by the joint dependence of the gas fraction in halos 
on the bulk vel ocity and the local overdensity in e ach pa tch. 

Following iTseliakhovich. Barkana fc Hiratal (|2010h we 
find the fraction of the baryon density contained in halos 
with mass larger than the minimum cooling mass M coo i 



/gas(> M coo l) 



M dn f g (M) 
Po dM f b 



dM 



(4) 



where po is the mean matter density today, dn/dM is the 
comoving abundance of halos of mass M, ft = f2i,/fi m is 
the mean cosmic baryon fraction and f g {M) is the fraction 
of the total halo mass which is in the form of gas. The 



1 Since u bc decays as 1 + 2 throughout the universe, a patch that 
has the root-mean-square value of v^, c at one redshift will have the 
root-mean-square value of the relative velocity at every redshift, 
and in particular v^ c = 30 km sec -1 at recombination. 




30 40 
Redshift 

Figure 2. The global mean gas fraction in star-forming ha- 
los (solid curves) and in minihalos, i.e., halos below the cooling 
threshold (dashed curves). The results, based on our optimal fit 
(cq. [3j are shown after averaging over the distribution of relative 
velocity (thick curves), or in the case of no relative motion, i.e., 
for V\, c (z) = (thin curves). 



gas fractions f g (M) depend on the filtering mass, which 
measures the scale at which the baryon fluctuations dif- 
fer substantially from those in the dark matter. In each 
patch, the filtering mass depends on the bulk velocity, and 
thus so do the gas fractions. Since the baryons contribute 
to the total power spectrum, the halo abundance dn/dM 
(which depends on fluctuations in the total matter den- 
sity) va ries as well with V hc- W e use the halo mass func- 
tion of ISheth fc Tormenl (ll999Tl . We refer the reader to 
ITseliakhovich. Barkana fc Hiratal (|2010l ) for the full details. 

We begin by recalculating some of the results of 
ITseliakhovich. Barkana fc Hiratal (|2010h . We show in Fig- 
ure [2] the redshift evolution of the globally averaged gas 
fraction in star-fo rming halos or in gas minihalos. Compare d 
with Figure 8 of ITseliakhovich. Barkana fc Hiratal (|2010h , 
our gas fractions are substantially lower, e.g., the gas frac- 
tion in halos above the minimum cooling mass is lower by 
a factor of ~ 3 at redshift z = 20. This is due to our 
higher A/ coo i and lower power spectrum normalization (see 
Section 6 for a full discussion of our differences with previ- 
ous papers). Note that the gas fraction in halos above the 
minimum cooling mass is proportional to the stellar mass 
density, assuming a fixed star formation efficiency (averaged 
over each 3 Mpc patch). 

In general, the importance of the relative velocities in- 
creases with redshift. Comparing the two categories of ha- 
los, we find that the relative suppression of the minihalos 
is larger than that of the star-forming halos at low red- 
shift; however, the relative suppression of the star-forming 
halos increases faster with redshift, and eventually it be- 
comes larger than that of the minihalos (beyond z ~ 50). 
At z = 20, the bulk velocities reduce the mean gas fraction 
in star-forming halos by a factor of 1.8 and that in minihalos 
by 3.1. 

Unlike previous studies, in our calculations the relative 
velocities produce three separate effects: suppression of the 
halo abundance (dn/dM), suppression of the gas content 
within each halo (f g (M)), and boosting of the minimum 
cooling mass (through V coo \(z)). In order to gain a better 
physical understanding, and for easier comparison with pre- 
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Figure 3. The ratio (compared to the v^, c = case) by which 
the bulk velocities change the global mean gas fraction in halos 
above the cooling mass (top panel) and in star-less minihalos 
(bottom panel). We consider four different cases: the full effect of 
the velocities (thick solid curves); the effect of v^ c in boosting the 
cooling mass only (dashed curves); the effect of i>b c in suppressing 
the halo abundance only (dotted curves); and the effect of v^ c in 
suppressing the gas fraction only (thin solid curves). 

vious papers, we investigate the relative importance of each 
effect in Figure[3] For the star-forming halos, the suppression 
of gas content is always the least significant effect (e.g., sup- 
pression by a factor of 1.13 on its own at z = 20), while the 
cooling mass boost is most important above z = 28.5 (factor 
of 1.26 on its own at z = 20), and the halo abundance cut is 
most important at lower redshifts (factor of 1.43 on its own 
at z = 20). For the minihalos, the boosting of the minimum 
cooling mass acts as a (small) positive effect, since it moves 
gas from the star- forming to the minihalo category (e.g., 
boost by a factor of 1.10 on its own at z = 20), while the 
other two effects are larger and comparable (e.g., at z — 20 
the suppression of gas content would give a reduction by a 
factor of 2.17 on its own, and the halo abundance cut would 
give a suppression factor of 1.74). 

4.2 Inhomogeneous gas fraction due to the 
dependence on the relative velocity 

The gas fractions shown in Figures 2 and 3 are globally 
averaged. However, in reality the universe is highly inhomo- 
geneous on small cosmological scales. We can divide it into 



patches that have various bulk velocities and densities. In 
this section we consider just the variation with velocity, i.e., 
averaged over all density fluctuations. In other words, we 
look at the contribution of velocity fluctuations to fluctua- 
tions in the gas fraction in halos. If we consider patches that 
are still small enough to have a coherent v\, c (e.g., cubes of 
3 comoving Mpc on a side), then the absolute value of the 
bulk velocity in each one follows a Maxwell-Boltzmann dis- 
tribution (eq. [TJ. 



Consider the contributions of patches of various veloc- 
ities to the total amount of star formation. At a given red- 
shift, the gas fraction in star-forming halos is lower in the 
patches with a high value of the relative velocity, because 
all three velocity effects (see the previous subsection) tend 
to reduce this gas fraction. On the other hand, patches with 
zero bulk velocity do not contribute much, simply because 
they are rare. As shown in the top panel of Figure [4] the 
most common bulk velocity is Ub c ~ 0.82<7„ bc , where Vbc and 
<r Vbc are both measured at the same redshift (recombination 
or any other z). If the stellar density were independent of 
the bulk velocity, then the contribution of regions of various 
velocities would be proportional to the velocity PDF. In- 
stead, the velocity suppression effect shifts the contribution 
to stellar density (assumed proportional to the gas fraction 
in star- forming halos) towards lower v^ c , with the relative 
change (compared to the Maxwell-Boltzmann distribution) 
increasing strongly with redshift. Thus, the biggest contri- 
bution to stellar density comes from v^ c = 0.67<7„ bc patches 
at 2 = 20, and from Vbc = 0.23<7t, bc patches at z = 60. We 
compare the contributions of the three separate effects of 
the velocity to the shift in the distribution of star formation 
(Figure[4] bottom panel). As in the top panel of Figure[3] we 
find that the suppression of halo gas content has the least 
significant effect on star-forming halos at z = 20 (typically, 
a ~ 10% effect on the distribution), while the other two ef- 
fects (halo abundance suppression and cooling mass boost) 
have a ~ 20 — 30% effect each. 



Thus, at the highest redshifts, the star formation is con- 
centrated in low- velocity regions which are rare, i.e., at the 
low-probability v^ c end of the Maxwell-Boltzmann distribu- 
tion function. The universe at these epochs is very inhomo- 
geneous, with a few bright regions filled with stars, while 
in all other regions the relative velocity is too high to allow 
significant star formation. As the universe expands, the rel- 
ative velocity decays, and in more and more patches across 
the universe the relative velocity drops enough to allow for 
star formation. As a result, the stellar distribution becomes 
increasingly homogeneous. To quantify the degree of inho- 
mogeneity caused by the dependence of stellar density on 
the bulk velocity, we plot the fraction of the volume of the 
universe (at lowest velocity, i.e., at highest stellar density) 
that contains 68% or 95% of the star-forming halos (Fig- 
ure [5}. The effect of volume concentration is mild at z = 20 
(68% of the stars are in 54% of the volume, and 95% in 
89% of the volume) , while it becomes very strong at z = 60 
(68% of stars in 4.6% of the volume, and 95% in 16% of the 
volume) . 
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Figure 4. Top panel: The relative contribution of regions with a 
given streaming velocity to the global gas fraction in halos above 
the cooling mass, i.e., df gas (> M coo i)/dv^, c normalized to an area 
of unity. The dependence is shown for z = 60 (solid curve) and 
z = 20 (dashed curve). We also show the Maxwell-Boltzmann 
distribution of the bulk velocity (dotted curve). The velocity is 
expressed in units of its root-mean-square value <x„ b . Bottom 
panel: The ratio at z = 20 between the quantity shown in the top 
panel (the relative contribution of regions with a given stream- 
ing velocity to the gas fraction in star-forming halos) and the 
Maxwell-Boltzmann distribution. If star formation were indepen- 
dent of bulk velocity, this ratio would equal unity. We consider 
this ratio for the same four cases as in Figure [3] the full veloc- 
ity effect (thick solid curve), the boost in the cooling mass only 
(dashed curve), the suppression of halo abundance only (dotted 
curve), and the suppression of the gas fraction only (thin solid 
curve) . 



4.3 Inhomogeneous gas fraction due to velocity 
and density fluctuations 

In order to quantify the full degree of inhomogeneity and 
concentration of star formation, we must include the effect 
of density fluctuations as well. In this section we thus con- 
sider the full PDF of the halo gas fraction within 3 Mpc 
patches, where the fluctuations result from a combination 
of the relative velocity distribution considered in the previ- 
ous section and density fluctuations. Specifically, the average 
density in a patch varies due to fluctuations on scales larger 
than its size. This average density follows a Gaussian distri- 
bution and is independent of the relative velocity within the 
same patch. 




Redshift 

Figure 5. The fractional volume of the universe that contains 
68% (dashed curve) or 95% (solid curve) of the star-forming halos 
as a function of redshift, where we consider just the contribution 
of velocity fluctuations to the inhomogeneity of star formation on 
3 Mpc scales. 



To find the modified halo mass function within a patch 
of a given overdensity 5r and bulk velocit y t^ c , we use the 
hybrid prescription (which combines the ISheth fc Tormenl 
( 1999) mass function with the extended P ress-Schechter 
model) introduced by Barkana fc Loebl d2004l) a nd general- 
ized by iTseliakhovich. Barkana fc Hiratal lj2010h to include 
Vbc- The dependence of the gas fraction in halos above the 
cooling mass on the two independent variables is illustrated 
in Figure [6] The dependence on both Sr and Vbc (each mea- 
sured in terms of its root-mean-square value) is stronger 
at higher redshifts. At a given redshift, the dependence on 
5r is stronger (i.e., the slope is higher) when v^ c is higher, 
since in this case the large halos (above the high cooling 
mass) are rarer and their abundance is more sensitive to 
the overdensity of the patch. If we consider the total range 
between and 2 a, we find that density and velocity fluctu- 
ations make comparable contributions to the star-formation 
fluctuations on the 3 Mpc scale. The relative importance 
of velocity increases with redshift and it will also increase 
if we consider larger scales. Even at z = 20 the velocity 
causes order unity fluctuations in the stellar density, and 
these fluctuations should be present at the large (100 Mpc) 
scales spanned by the velocity correlations. 

The resulting full PDF of the halo gas fraction is shown 
in Figure [7] (top panel), both for the star-forming halos, and 
the star-less gas minihalos. The main effect of the bulk veloc- 
ities is to shift the distributions towards lower gas fractions. 
At redshift 20, the effect is larger on the minihalos. In Fig- 
ure [7] (bottom panel) we show the fraction of the volume of 
the universe (at the high gas fraction end of the full PDF) 
that contains 68% or 95% of the stars, with and without the 
velocity effect. 

The volume concentration of star formation is a result 
of a complex interplay of the two sources of fluctuations. 
The global star formation is highest in the rare regions with 
both low bulk velocity and high overdensity, but more gen- 
erally, one of these can compensate for the other. The effect 
of Vhc on star-forming halos vanishes by z ~ 10, in agree- 
ment with our previous results, leaving just the effect of 
the local density. Even at somewhat higher redshifts (up to 
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Figure 6. The percentage of the gas fraction in star-forming 
halos at redshifts z = 20 (thick curves) and z = 40 (thin curves) 
as a function of the average overdensity 5r in the 3 Mpc patch 
(normalized by its root-mean-square value ctr), for various values 
of the relative velocity: no relative motion (dashed), v^, c = cr Vbc 
(solid) and v^ c = 2cr Vbc (dotted). 



z ~ 35), the concentrating effect of the velocities on their 
own (Figure [SJ) remains weaker than that of the densities 
alone (no- velocity case in Figure U), so at these redshifts the 
full case is dominated by the densities, and the concentrat- 
ing effect of density is enhanced by including the velocities 
(which steepen the dependence on density: Figure[6jl . At red- 
shifts above ~ 35, velocities dominate, and then including 
the density fluctuations (compared to averaging over them 
at each velocity) actually reduces the concentration since 
it allows low-velocity regions to contribute relatively more 
volume with high gas fractions (due to the steeper density 
dependence at high bulk velocity). 

Specifically, at z = 20, density fluctuations alone (i.e., 
setting i>bc = 0) would concentrate 68% of the stars into 
39% of the volume and 95% into 81% of the volume. The 
addition of the bulk velocity provides a mildly increased con- 
centration into 35% and 77% of the volume, respectively. 
At redshift 60 the results are that 68% of the stars are in 
11% of the volume and 95% in 45% (which is higher than 
in Figure [5]), compared to 14% and 52% of the volume, re- 
spectively, at zero bulk velocity. The effect of the velocities 
should be more clearly apparent on scales larger than our 
3 Mpc pixels, i.e., in addition to the small additional con- 
centration that they cause (as seen in Figure [7]), their effect 
is to redistribute the star-forming regions to produce larger 
coherent regions of either high star formation or low star 
formation (voids). 

We note that the assumption that the local overden- 
sity on large scales 5r and the streaming velocity «bc are 
statistically independent is not perfectly accurate. A patch 
with a high local overdensity has expanded less than other 
patches, so that the peculiar velocity Vbc has not declined 
as much compared to the expansion. Indeed, we expect that 
Vbc fb c (l + 5r/3). However, we have found that this cor- 
rection makes only a small difference to the PDF (up to a 
4% relative error at z = 60, and less at lower redshifts). 
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Figure 7. Top panel: The full probability distribution function 
(PDF) of the gas fraction at redshift z = 20. We show the PDF 
of the gas fraction in halos above the cooling mass (solid curves) 
and the PDF of the gas fraction in star-less minihalos (dashed 
curves). We consider two cases: randomly distributed u bc and 8r 
(thick curves), and v^c = but random 8r (thin curves). Bottom 
panel: The fractional volume of the universe that contains 68% 
(dashed curves) and 95% (solid curves) of the star forming halos, 
where we consider the full PDF in 3 Mpc patches. In each case 
we consider including the relative motion (thick curves) or not 
(^bc = 0, thin curves). 



5 THE FIRST STAR 

In the previous sections we have discussed the conditions 
needed to initiate star formation. The main condition is that 
the halo mass must be large enough to allow molecular cool- 
ing. Given a large enough initial density fluctuation, a halo 
with a sufficiently large mass will form relatively early. The 
very first stars depend on extremely rare fluctuations, hence 
we need to average over the volume of the observable uni- 
verse, (14 Gpc) 3 , in order to have the full statistical range 
needed to accurately estimate the formation time of the first 
star. 

Due to computational limitations, numerical simula- 
tions can form sta rs only in a very lim ited cosmological con- 
text. For instance. iGreif et alj (1201 ll) studied s tar formation 
in a (500 kpc) 3 volume and lStacv et alj l|2010l ) were limited 
to (100 h _1 kpc) 3 . In a small volume the chance of getting 
a rare high density fluctuation is quite small. Therefore the 
formation redshift of the first stars in simulations is greatly 
underestimated, with most simulations forming their first 
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star below redshift 30. The highes t redshift where a star has 
formed in a simulation is z = 47 dReedet al.ll2005l ). 

iNaoz. Noter fc Barkanal (|2006l ) first applied these sta- 
tistical considerations in order to predict the redshift of the 
first observable star (i.e., in our past light cone) analyti- 
cally. They estimated the redshift of the first star to be 
z = 65, using the 3-year WMAP set of cosmological parame- 
ters (|Spergel et alj|2007h and assuming a minimum circular 
velocity for cooling of V CO oi = 4.5 km sec -1 . In this section 
we generalize their method in order to account for the bulk 
velocities and estimate their impact on the epoch of the first 
star formation. This problem is particularly relevant since 
the effect of the relative velocity on star formation increases 
with redshift, and is thus at its maximum when we con- 
sider the very first star. We also study the sensitivity of the 
first-star reds hift to various uncertainti e s. 

Following INaoz. Noter fe Barkanal l|2006h we calculate 
the mean expected number (N(> z)) of star-forming halos 
that formed at redshift z or higher, but where the halo abun- 
dance is now averaged over the bulk velocity distribution at 
each redshift. This number is the ensemble-averaged number 
of stars, but we have only one universe to observe. Hence, 
we expect Poisson fluctuations in the actual observed num- 
bers. The probability of finding at least one star is then 
1 — exp [— {N(> z))], and (minus) the redshift derivative 
of this gives the probability distribution p*(z), where the 
probability of finding the first star between z and z ± dz is 
p,(z)dz. 

As shown in Figure [8] (top panel) , we find that in the 
absence of the bulk velocities, the first star would be most 
likely to form at z = 69.9, with a median z = 70.3 (cor- 
responding t o t = 30 Myr after the Big B ang). The dif- 
ference with INaoz. Noter fc Barkanal (|2006h is due to the 
changes in the cosmological parameters between WMAP3 
and WMAP7, specifically the increased power on the rele- 
vant scales (since the increased spectral index has a larger 
effect than the reduced as), and the decreased cooling mass 
in the Vhc = case compar ed to the value assumed by 
INaoz. Noter fc Barkanal <|2006h . 

The relative velocity effect delays star formation, where 
for the very first star we find a delay of Az = 5.3. The first 
star is now most likely to form at z — 64.6, with a median 
z = 65.0 (t = 34 Myr) that has a 1 - a (68%) confidence 
range z = 63.9 — 66.5 due to the Poisson fluctuations. In 
addition, the redshift of the first star is uncertain due to 
the current errors in the cosmological parameters and the 
uncertanity in the cooling mass. Regarding the cosmologi- 
cal parameters, the redshift of the first star is sensitive to 
the amount of power o n the scale of the fir st halos. The 
uncertainly of WMAP7 (jKomatsu et al.ll2010l ) in the ampli- 
tude of the primordial fluctuations (parameterized by as) is 
A(J8 = ±0.024, which implies (for our optimal fit) an un- 
certainty of Az — ±2.2 in the median redshift of the first 
star. The larger is as, the earlier will the first star form. 
More generally, we include the current correlated errors in 
the full suite of standard cosmological parameters, and find 
a resulting Az = ±5.1. 

In order to estimate the impact of the current uncer- 
tainty in the effect of the bulk velocity on the minimum 
cooling mass, we estimate the redshift of the first star for 
each of the fits discussed in Section 3. We find (Figure [8] 
bottom panel) that the range of the SPH simulations is a 
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Figure 8. Top panel: The impact of the relative velocity on the 
redshift of the very first observable star. We plot the probabil- 
ity density of seeing the first star at a given redshift, including 
the effect of relative velocity for our optimal fit (solid curve), or 
without the effect of the velocity (i.e., for the same fit but with 
v bc = 0, dotted curve). The formation of the first star is delayed 
by Az = 5.3 due to the relative velocity effect. We mark the me- 
dian redshift of the first star for each distribution (•). Bottom 
panel: The probability density of the redshift of the first star 
calculated for each of the fits of Figure [T] The median redshifts 
of the first star (from left to right ) are: z = 63.5 ( "fit" to the 
AMR simulations), z = 64.3 (fit to lGreif et ah! feoill) '). * = 65.0 
(the optimal fit t o the SPH simulations) and z = 65.8 (fit to 
IStacv et al.1 teulOh l. 



Az — 1.5, and the discrepancy between the AMR and SPH 
simulations is comparable. Thus, we conclude that the delay 
due to the bulk motion is substantial, but there are still sig- 
nificant uncertainties in it. In summary, we find the median 
redshift of the first star in our observable universe to be 



z = 65. 0_ 1 'i (Poisson) l r (simulations) ± 5.1 (cosmology) . 

(5) 

Thus, current uncertainties in the values of the cosmological 
parameters dominate over the differences in the simulations 
and the irreducible Poisson fluctuations. 



6 DISCUSSION 

We have studied the impact of the relative motion between 
the gas and the dark matter on the formation of the first 
stars. We included a new effect found in recent small-scale 
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hydrodynamic simulations. In particular, we fit their results 
to a physically-motivated ansatz that expresses the mini- 
mum circular velocity of gas-cooling halos as a simple func- 
tion of the local bulk velocity when the halo forms. This 
result implies that in contrast to previous expectations, the 
minimum mass of star-forming halos does not decrease with 
redshift, except in regions with very low values of the bulk 
velocity. 

This result implies that the relative velocities produce 
three separate effects: suppression of the halo abundance, 
suppression of the gas content within each halo, and boost- 
ing of the minimum halo mass required for cooling. Quanti- 
tatively, we found that the halo abundance cut has a large 
effect on the two categories of halos (star-forming halos and 
star-less minihalos), while the cooling mass boost primarily 
affects star-forming halos and the suppression of gas content 
primarily affects the minihalos. In total, at z = 20 the bulk 
velocities reduce the mean gas fraction in star-forming halos 
by a factor of 1.8 and that in minihalos by 3.1. Thus, even 
at 2 = 20 the velocity causes order unity fluctuations in the 
stellar density, and these fluctuations should be present at 
the large (100 Mpc) scales spanned by the velocity correla- 
tions. 

The velocity dependence of the gas fraction tends to 
concentrate the global star formation into regions of low 
bulk velocity. In particular, at z — 20, 68% of the stars are 
in the 54% of the volume with the lowest velocity, and 95% 
are in 89% of the volume. Adding in the effect of density 
fluctuations tends to concentrate the global star formation 
into regions of both low bulk velocity and high overdensity. 
As a result, at z = 20, 68% of the stars form within 35% of 
the volume and 95% in 77% of the volume. This concentra- 
tion effect becomes much stronger at higher redshifts. 

The formation of the very first star is delayed by Az ~ 5 
due to the bulk velocities. Given the updated cosmologi- 
cal and astrophysical parameters, the first star is now most 
likely to form at z = 64.6, with a median z — 65.0 (corre- 
sponding to a cosmic age of t = 34 Myr). Due to the com- 
bination of density and velocity fluctuations, the formation 
of stars begins at different times in different regions. This 
leads to a very inhomogeneous early universe. Although by 
redshift 20 most of the universe is populated, the age of the 
oldest stars in each region is significantly different. 

To make the novelty of our work clear, we now 
make a full comparison of the ingredients of our calcula- 
tions with those in the previ ous literature. We start with 
iTseliakhovich fc Hiratal l|201Cf ). who discovered that the rel- 
ative velocity effect is important. They only calculated the 
impact on the halo abundance, but this was sufficient for 
them to deduce the implication of large-scale fluctuations. 
However, their calculations had a number of simplifying as- 
sumptions: they calculated the baryon perturbations under 
the approximation of a uniform sound speed, and used the 
old Press-Schechter halo mass function. 

iDalal. Pen fc Seliakl |2010h were the first to point out 
the effect of the relative velocity on suppressing the gas con- 
tent of halos. However, they made a number of simplifying 
approximations that we have relaxed here. These include: 

(i) We have calculated the filtering mass (Mf) from 
linear theory, while they took the effective value found in 
simulations in the standard (no relative velocity) case, and 
then multiplied it by a simple Wbc-dependent ansatz. 



(ii) We have allowed for a smooth transition between 
gas-rich halos atM> Mf and gas-poor halos at M < Mf 
as is suggested by simulations, rather than applying a step- 
function cutoff. 

(iii) We have simultaneously included the dependence 
of the gas fraction in halos on the large-scale matter over- 
density 8r and relative velocity Vbc- This combines both 
the "traditio nal" biasing model (which i ncludes 5r but not 
«bc) an d the IDalal. Pen fc Seliakl (|2010t ) treatment (which 
includes Vbc but not Sr). We found that both effects are 
important (compare Sections 4.2 and 4.3). 

(iv) We incl uded the effect of v b c on the halo 
mass function (ITseliakhovich fc Hiratal l201(t) . which 
IDalal. Pen fc Seliakl ||201Ch did not. ~ 

(v) Most importantly, we incorporated a cooling crite- 
rion for star formation, rather than scaling by the total gas 
content in halos. The vast majority of the gas is in minihalos 
that cannot cool, and because of their low circular velocities 
their ability to collect baryons is much more affected by «b c 
than the star-forming halos. This suggests that the effect of 
relative velocities on early sta r form ation might be less than 
found bv lDalal. Pen fc Seliakl (|2O10Tl . However, we find that 
the inclusion of the other effects (mass function and cool- 
ing threshold, in addition to baryon fraction) does restore 
the expectation for order unity fluctuations, with exciting 
implications for observational 21-cm cosmology. 

In part of this paper we closely followed 
ITseliakhovich. Barkana fc Hiratal (|2010l 'l. However, we 
fixed two inaccuracies in their power spectrum (in the nor- 
malization and the spectral slope) that gave substantially 
too much power on small scales. Then, our main goals 
were to include the new effect on the cooling mass based 
on simulations, to extend the calculations to the highest 
redshifts of star formation, and to quantify the degree of 
concentration of star-forming halos. With there now being 
three separate effects of the bulk velocity, we also carefully 
studied the relative importance of these various effects. 
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